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Certain aspects of the recently proposed antisymmetrised a particle product state wave function, 
or THSR a cluster wave function, for the description of the ground state in **Be, the Hoyle state in 
^^C, and analogous states in heavier nuclei, are elaborated in detail. For instance, the influence of 
antisymmetrisation in the Hoyle state on the bosonic character of the a particles is studied carefully. 
It is shown to be weak. Bosonic aspects in Hoyle and similar states in other self-conjugate nuclei are, 
therefore, predominant. Other issues are the de Broglie wave length of a particles in the Hoyle state 
which is shown to be much larger than the inter-alpha distance. It is pointed out that the bosonic 
features of low density a gas states have measurable consequences, one of which, that is enhanced 
multi-alpha decay properties, likely already have been detected. Consistent with experiment, the 
width of the proposed analogue to the Hoyle state in ^^O at the excitation energy of Ex = 15.1 
MeV is estimated to be very small (34 keV), lending credit to the existence of heavier Hoyle-like 
states. The intrinsic single boson density matrix of a self-bound Bose system can, under physically 
desirable boundary conditions, be defined unambiguously. One eigenvalue then separates out, being 
close to the number of a's in the system. Differences between Brink and THSR a cluster wave 
functions are worked out. No cluster model of the Brink type can describe the Hoyle state with a 
single configuration. On the contrary, many superpositions of the Brink type are necessary, implying 
delocalisation towards an a product state. It is shown that single a particle orbits in condensates 
of different nuclei are almost the same. It is thus argued that a particle (quartet) antisymmetrised 
product states of the THSR type are a very promising novel and useful concept in nuclear physics. 



I. INTRODUCTION 

Recently it has been pointed out that certain states 
in self conjugate nuclei around the alpha-particle disinte- 
gration threshold can be described as product states of a 
particles, all in the lowest QS state. Considerable theo- 
retical and experimental activity has developed since this 
idea was first put forward in 2001 In this paper we 
want to further dwell on the basic foundations and pre- 
dictions of this concept, since the usefulness of the latter 
has recently been questioned 0]. 

In Refs. d, 0, Q it was pointed out that in homoge- 
neous nuclear matter a-particle condensation is a pos- 
sible nuclear phase. Therefore, the above mentioned a- 
particle product states in finite nuclei have been proposed 
to be related to boson condensation of a particles in infi- 
nite matter The infinite matter study used a four 
particle (quartctting) generalisation of the well known 
Thouless criterion for the onset of pairing as a function 
of density and temperature. The particular finding in the 
four nucleon case was that a-particle condensation only 
can occur at very low densities where the quartets do not 
overlap appreciably. This is contrary to the pairing case 
where, in weak coupling situations, the Cooper pairs also 



may strongly mix. It is interesting to note that the low 
density condition for quartetting was in the meanwhile 
confirmed in Ref. |^ with a theoretical study in cold atom 
physics. 

Concepts developed for infinite nuclear matter are of 
value also to interpret properties in finite nuclei and to 
construct useful approximations. As examples, we re- 
fer to pairing, two and more body correlations, and one 
body occupation numbers. Pairing is believed to occur 
in neutron stars which are considered as infinite neutron 
matter becoming superfluid below a critical temperature. 
Pairing also is a useful concept in many finite nuclei, in 
spite of the fact that nuclei are not macroscopic objects. 
Therefore, in reality, they are only in a fluctuating state 
and we have to project, e.g., a BCS state on a definite 
number of nucleons. In spite of the finiteness of nuclei, 
the BCS state remains a useful approximation for the 
quantum state. For example, the strong reduction of 
measured moments of inertia of such nuclei compared 
with the classical values are explained as a consequence 
of superfluidity 0] ■ 

Alpha-particle condensed states may also be of rele- 
vance in finite nuclei. As already pointed out above, in 
Ref. [l|, see also Q, we interpret the Hoyle state, i.e. 
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the Oj-state at = 7.65 MeV in ^^C, as a product 
state of three a's and predict that Hoyle-Hke states very 
Ukely also exist in low density states in heavier na nu- 
clei, close to the na disintegration threshold. Examples 
for ^^O and ^"Ne have been presented in Refs. [H, @, em- 
ploying the so-called THSR wave function, proposed by 
Tohsaki, Horiuchi, Schuck, and Ropke. It is analogous 
to the (number projected) BCS wave function, replac- 
ing, however. Cooper pairs by a-particles (quartets). In 
addition, we showed that pure product states of a parti- 
cles in the threshold states are realised to about 70 per- 
cent [0, 113 • We, thus, define a state of condensed 
na's, if in a nuclear state the latter forms in good approx- 
imation a bosonic product state, all bosons occupying the 
lowest quantum state of the corresponding bosonic mean 
field potential. 

In the present paper we will further investigate the 
concepts and consequences of the THSR wave function. 
We will carefully study the effect of antisymmetrisation 
on the bosonic character of the a particles in the Hoyle 
state. It will be found that, compared to the ground 
state, its influence is very weak, however, not negligible 
and in a sense necessary for the a particle gas state of 
low density to be formed and stabilised. We also focus on 
measurable properties that are adequately described in 
this approximation. Considering special results we show 
the usefulness of this novel type of a cluster wave func- 
tion. 

The paper is organised as follows. In Sec. II, we elab- 
orate on the difference between Brink and THSR type 
wave functions and demonstrate this first for the most 
simple case of ^Be. Sect. HI is dedicated to a study of the 
antisymmetrisation effects between the a's in the Hoyle 
state and in Sec. IV the de Broglie wave length of a's 
in the Hoyle state is studied. In Sec. V we discuss some 
measurable consequences of the bosonic character of a's 
in compound states. Sec. VI considers decay properties 
of Hoyle-like states, and Sec. VII reveals the similarity 
of a wave functions in Hoyle-like states in nuclei with 
different numbers of a particles. Finally in Sec. VIII, we 
conclude. 



II. DELOCALISED ALPHA-PARTICLE 
CONDENSATE VS LOCALISED BRINK TYPE 
WAVE FUNCTIONS 

In Ref . 01 , it is claimed that the proposed THSR wave 
function for the description of loosely bound a particle 
states is an approximation to existing nuclear a-cluster 
states. The authors essentially have in mind localised 
cluster states of the Brink type In this section we 
will show in detail, presenting new investigations, that 
the situation is just the contrary: states which have well 
born out na-cluster structures like, e.g. the ground and 
low lying states of ^Be and the Hoyle state in ^^C are 
much more adequately described by THSR than by Brink 
wave functions. 



Wave functions for self conjugate light nuclei which in- 
corporate a-cluster substructures have by now been in 
use in nuclear physics since about half a century p^ . 
Two na nuclei have been on the forefront of the investi- 
gations: ^Be and ^^C, with two and three alpha particles, 
respectively. The starting point always has been practi- 
cally the same, that is for the individual alpha parti- 
cles an intrinsic translationally invariant mean field wave 
function of Gaussian type, representing the free space 
alpha-particle wave function is taken, whereas the cen- 
ter of mass (c.o.m.) motion is either determined from a 
full variational principle or limited parametrised ansatze 
reflecting certain underlying ideas of the alpha-particle 
motion have been assumed. Let us, therefore, write for 
an na nucleus the typical following alpha-cluster wave 
function 

X(?!'Qi(ri,i, • • • ,ri4) • • •0a„(r„,i, • • • ,r„,4)], (1) 

with A the antisymmetriser, Ri the c.o.m. coordinates 
of the a-particle, and 
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the normalised intrinsic wave function of the a-particle 
with (OS')"' shell-model configuration. With b a varia- 
tional parameter, it is well known that Eq. ([2]) can de- 
scribe the free a-particle quite well. The a-particle wave 
functions in Eq. ([1]) arc fixed to their free space form. 
The wave function x for the c.o.m. motion of the a's 
with Ri = jivi^i + ri^2 + '"i,3 + ^1,4) is, of course, also 
chosen translationally invariant, that is it depends only 
on the relative coordinates Rij = Ri — Rj or on the cor- 
responding Jacobi coordinates. The spin-isospin part in 
Eq. ([T]) is not written out but supposed to be of scalar- 
isoscalar form. We will not mention it hitherto. As al- 
ready pointed out, the c.o.m. wave function is either 
determined from a full variational calculation, minimis- 
ing the energy, or one adopts a restricted variational 
ansatz. A famous example is the so-called Brink wave 
function which places the a-particlcs at certain positions 
in space [l^. For example in the case of *Be, this is 

x|""''(-Ri,H2) ocexp [-(^Ri-R/2-iR2+R/2)y /b^ 

with an obvious generalisation to the case of na-particles. 
Usually in Eq. ^ b has the same value as for the free 
space a-particles and then Eq. ([3|) implies that the two 
a's arc placed a relative distance R apart from one an- 
other. Though this kind of geometrical, crystal-like view 
of the cluster structure works well for many cases, for 
instance, parity-violating ^^C-|-a, ^^0-|-a, and '^''Ca-l-a 
structures in ^^O, ^°Ne and ''^Ti, respectively [isl . [l6llT7t . 
and also when additionally neutrons are involved [l8[ , it 
is on the contrary known since several decades that this 
picture fails for the description of the famous Hoyle state, 
i.e. the Oj state in ^^C. 
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Since such basic results of cluster physics may not be 
common knowledge, we here want to present a study with 
Brink-type cluster wave functions and also compare it 
with another variational ansatz for x with the diamet- 
rically opposite point of view of completely delocalised 
a-particles, namely the THSR wave function cited above 
There, the following form for x is adopted 

x™^^(fii, i22, • • • , Rn) = xo(7^l)xo(7^2) • • • xo(7^„), 

(4) 

with Tti = Ri - Xg where Xq = (-Ri H h Rn)/n is 

the total c.o.m. coordinate, and 



Xo(7^) = exp[-2(-RVs')], 



(5) 



that is a Gaussian with a large width parameter B which 
is of the nucleus' dimension. The product of n identical 
OS' wave functions reflects the boson condensate charac- 
ter, discussed in the introduction. This feature is realised 
as long as the action of the antisymmetriser in Eq. ^ is 
sufficiently weak. It also is useful to notice that with @ 
and ([T]), we can write for Eq. ([1]) 



• • • ,r„,4|THSR) = A[lpa^i^a 



with 



and 



|THSR) = |THSR(B)) = A\B) 



(^1,1, • • • , rn,i\B) = tpat^a2 ' ' ' 4'a„ 



(6) 



(7) 



(8) 



where ipa^ = xo{'^)4'ai and the definitions Eqs. ([7]) and 
© will be useful below. Equations ©, ^ and (g]) high- 
light the analogy of the THSR wave function with the 
number projected BCS wave function for pairing 

(ri_i, • • • ,r„,2|BCS) = A[(j)pi,h-{ri^i,ri^2) 

■ ■ ■(f>paiiirn,i,r„,2)], (9) 

where (/)pair(''i,i, '"1.2) is the Cooper pair wave function. 

This type of condensate wave function has known, in 
the meantime, considerable success, notably with an ac- 
curate description of the Hoyle state, proposing it as the 
first of a series of excited states in na nuclei with a- 
particle product character. Usually one employs Jacobi 
coordinates and then the THSR ansatz for x is given 
by: 



THSR 
Ana 



exp 



' P + 2/32 



(10) 



with Hi ~ + A slight generalisation of (|10p is pos- 
sible, taking into account nuclear deformation, see below. 

Before discussing the Hoyle state, we want to study ^Be 
in some detail because even this nucleus which is known 
to have intrinsically a two-alpha dumbbell structure can 
very well be described in the laboratory frame with the 



delocalised THSR wave function. Let us repeat Eq. ([T|) 
for this particular case 



^2a = A[x{Rl2)(f>ai(l)a2], 



(11) 



with R12 = Ri — R2- Note that (|lip is a fully anti- 
symmetric and translationally invariant wave function in 
8—1 = 7 coordinates. Minimising for a given Hamilto- 
nian with N-N forces and Coulomb force [lj| , the ground 
state energy with respect to x, leads straightforwardly to 
a Schrodingcr type of equation for Xj corresponding to 
the Resonating Group Method (RGM) [li,[20|. 



K (f>ajH~ E\A[xir)(f>a., 0a J) = 0. 



(12) 



With the usual definition of RGM, this equation is trans- 
formed into a standard Schrodingcr equation for the wave 
function '^2a{''') of the relative motion of the two a- 
particles 



dr'/i(r,r')*2a(r') -^*2a(r) 
«'2a(r) = 



h{r,r') = 

where 

n{r, r') 
h{r,r') 



(13) 

dr'n^'^{r,r')x{r'), (14) 

dridr'^n~^l^{r,ri)h{ri,r[)n-^/'^{r[,r'), 

(15) 

{S{R,2^r)<ll\A[5{Rr2-r')^l]), 

{5{R,2 - r)cl,l\H\A[S{R^2 - r')0^]). (16) 



In Eq. (|16p . H is the microscopic Hamiltonian under 
consideration and 4>ai4>a2 is abbreviated by (j)^. 



Full (RGM) 
Brink (iJ=3.45) 




15.0 



r[fm] 



FIG. 1: Radial parts of wave functions '^2a{r) for the full 
RGM solution (full line) and single "Brink" component with 
R = 3.45 fm (dotted line). The Volkov force No.l is taken 
with Majorana parameter value M = 0.56. 



As mentioned, Eq. (fT^ and ([T5)) have been solved with 
very high accuracy since 50 years with excellent results 
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for all low energy properties of ^Be [1^ . The radial part 
of the wave function ^'2q;(t') is shown in Fig. [1] by the full 
line. We see that there exist two nodes, an effect which 
stems from the Pauli principle. We now will discuss two 
approximate forms for x(r) which are based on the two 
diametrically opposite views of the nature of ®Be already 
mentioned above: the THSR wave function and the Brink 
cluster wave function. Let us start with the latter. We 
have from Eq. ([3]) 



P 



j=o 



exp 



R) 



(17) 



with p-^=^ the projection operator on J = 0. In Eq. p7|) 
i? is a parameter which allows to place the two a's a 
distance R apart, and r is the relative coordinate be- 
tween the two a's, i.e. r — Ri — R2. This ansatz seems 
reasonable, since the microscopic calculation of Ref. [2l| 
indeed indicates that the two a's are about 4 fm apart. 
Obviously, the parameter R can be varied to find the 
optimal position of the a-particles. The result of such a 
procedure is shown on Fig. [T] with the broken line, taking 
the optimal value R = 3.45 fm (6 is kept fixed at its free 
space value, h = 1.36 fm). Qualitatively such a "Brink" 
wave function follows the full variational solution (full 
line). However, in the outer part, for instance in the ex- 
ponentially decaying tail quite strong differences appear. 
The squared overlap with the exact solution is 0.722. Of 
course, such Brink wave functions also can serve as a basis 
and it is interesting to study the convergence properties. 
We, therefore, write for the ^Be wave function appearing 
in Eq. (ini) 



-^[X(^)0ai(/'o 



fAir,R('\b) (18) 



and 



(19) 



where the i?*-'-* indicate the various positions of the a- 
particles and fi are the expansion coefficients. The con- 
vergence of the squared overlap with the exact solution 
is studied where we take for the positions R,^^^ = 1 fm, 
= 2 fm, • ■ •, = n fm. We start with n = 5. In 
Fig. m the convergence rate is shown as a function of n 
of the squared overlap with the exact solution and same 
for the energy. The point of n = 1 is with the optimised 
single Brink wave function (i?(i) = 3.45 fm). We see that 
the convergence is not extremely fast but for n = 20 the 
squared overlap with the full RGM solution amounts to 
0.9999. Also energy is converged to within 10~^. In Fig. [3] 
we show the convergence of the wave function r'^2aif)- 
In the insert we see that there is still a slight change in 
the far tail going from n = 25 to n = 30. 

Let us now investigate the THSR ansatz for x(r). 
There it is assumed from the beginning that the a's are 
delocalised and a single Gaussian ' centered at the 
origin with, however, a large width = + with (3 
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FIG. 2: Binding energy given by the wave function <l>2(i and 
the squared overlap between '^2a and $2q. For n = 1, a single 
Brink wave function with optimised R = 3.45 fm is adopted. 
See text for definition of $20 and ^2a. The Volkov force No.l 
is taken with Majorana parameter value M = 0.56. 



a variational parameter, is taken. Very much improved 
results over the single component Brink wave function 
are obtained. With (3 ~ 3.24 fm the squared overlap 
becomes 97.24%. However, practically 100 % accuracy, 
compared with the exact solution, can be achieved start- 
ing with a slightly improved ansatz, i.e. with an axially 
symmetric deformed Gaussian which is then projected 
on the ground-state spin J = (projections on J = 2, 4 
yield the rotational band of ^Be) |22| . 



THSR 
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P^=° exp 
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exp( 



-rV^i) 



52 

Erf(i 



2/31 



62 + 2/32 

.^2)1/2 



(20) 



and Erf(a;) the 



with Bf = 62 _|_ 2/32 and r\ = rl ^ 
error function. The second line of Eq. (j20|l is obtained 
from a simple calculation. 

Such an intrinsically deformed ansatz is, of course, 
physically motivated by the observation of the rotational 
spectrum of ^Be indicating a large value of the corre- 
sponding moment of inertia. The minimisation of the 
energy yields (3^ = Px = Py = 1-78 fm and /3z = 7.85 
fm. With these numbers, the squared overlap between 
the exact '^2a and 5*™^^ is with 0.9999 extremely pre- 
cise [i^l ■ In Fig. [3] we also show that the THSR wave 
function agrees almost completely even far out in the tail 
with the "exact" solution with 30 "Brink" components. 
On the scale of the figure exact and THSR wave func- 
tions cannot be distinguished. Let us also mention that 
beyond r ~ 3.5-4.0 fm t^™sr xj>^^ become practically 
identical in shape, except for a difference of normalisa- 
tion, meaning that approximately from the maximum on, 
the a-particles are not influenced any longer by antisym- 
mctrisation and behave as pure bosons. 
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FIG. 3: Comparison of THSR wave function with single com- 
ponent "Brink" wave function {n = 1). The convergence rate 
with the superposition of several (n) "Brink" wave functions 
is also shown (see text for more details). 
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FIG. 4: Probability distributions |r*2c(r)|^ for the THSR 
and single component "Brink" (denoted as n = 1) wave func- 
tions and for the full RGM solution. On the scale of the fig- 
ure the two curves of THSR and RGM are indistinguishable. 
They are normalised as J dr\r'^ 2a{r)'^ = 1. 



heavily overlap, because of the strong action of the Pauli 
principle, the relative wave function has two nodes and 
small amplitude, as shown in Fig. [3) Contrary to the 
outer part of the wave function determined dynamically, 
the behavior of the relative wave function in this strongly 
overlapping region is determined kinematically, solely re- 
flecting the r-dependence of the norm kernel in Eq. (jl6|) . 
This is clearly seen from the fact that both THSR and 
Brink wave functions have very nearly the same behavior 
in this region. It is also instructive to show the probabil- 
ity |r\['2a('')P which is presented in Fig. U) We see that 
the latter is practically zero for r < 2-3 fm, reminiscent of 
the excluded volume picture which is sometimes ado pted 
when the alpha-particles are treated as pure bosons [2J] . 
Let us repeat: the alpha's in *Be move practically as 
pure bosons in a relative 05 state of very low frequency 
as long as they do not come into one another's way. that 
is as long as they do not overlap. One should stress that 
this picture holds after projection on good total momen- 
tum and good spin, that is in the laboratory frame. It is 
equally true, as already mentioned, that in the intrinsic 
frame *Be can be described as a strongly deformed two 
alpha structure, sec ansatz (|20p . reminiscent of a dumb- 
bell. 
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FIG. 5: 



Structure of the 



state in C shown by the over- 



As we saw above, the single component, two parame- 
ter THSR ansatz, Eq. (|^D|) . for the relative wave function 
of two alpha's seems to grasp the physical situation ex- 
tremely well. The most important part of this wave func- 
tion is the outer one beyond some 3 fm. There, the two 
alpha's are in an S wave of essentially Gaussian shape. 
The corresponding harmonic oscillator frequency is esti- 
mated to fiio ^ 2 MeV. Therefore, as long as the two 
alpha's do not overlap strongly, they swing in a very low 
frequency harmonic oscillator mode in a wide and de- 
localiscd fashion, reminiscent of a weakly bound dimer 
state. Inside the region r < 2-3 fm where the two alpha's 



lap between the standard Brink cluster wave function of the 
isosceles configuration and the exact wave function. Figure 
is adopted from Ref. [25j . 



For the Hoylc state, it is known since long that the 
situation is qualitatively similar with, however, a much 
reduced action of the Pauli principle (see discussion be- 
low). In the work by Uegaki et al. [H, a contour 
plot of the overlap between a Brink type wave function 
and the full RGM solution for x of Eq. ([1]) is shown. 
We reproduce this figure in Fig. (5] It is seen that the 
overlap between the standard cluster wave function and 
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the exact solution is quite poor. In the best case the 
squared overlap reaches only about 50%. The authors 
of that article which dates three decades back, in view 
of their finding in Fig. [5l make the following significant 
statement to characterise the situation: ". . . the J state 
has a distinct clustering and has no definite spatial con- 
figuration. In other words, C is dissociated in the J 
state into weakly interacting three a-clusters, which can 
be considered as a three boson system The situa- 

tion also is highlighted in a recent work by Chernykh et 
al. [1^ where about 55 components of the "Brink" type 
wave functions are needed to reproduce the full RGM so- 
lution for the Hoyle state accurately, that is considerably 
more than in the case of ®Be. 

In what concerns the THSR wave function for the de- 
scription of the Hoyle state, the situation is slightly more 
complicated by the fact that the loosely bound three a 
configuration is now no longer the ground state but the 
state at 7.65 McV excitation energy We, therefore, 
have to minimise the energy with the THSR wave func- 
tion 



exp 



2 

1=1 



-52 



2/31 



52 + 2pl 



(21) 



where the 2 ^^'^ the two Jacobi coordinates, and /ii = 
1/2, ^2 = 2/3, under the condition that the 3a wave 
function $30, is orthogonal to the ground state. $30. to 
be used for the minimisation of the width parameters can 
schematically be written as 



A 



(22) 



where pj^ '^ -* keeps (p2)) orthogonal to the ground state 
configuration, i.e. pj^ '^ -' = l — |0f )(0]^|. The wave func- 
tion thus obtained has 99.3% squared overlap p7j with 
the full RGM solution of Kamimura et al. [2a| • The cor- 
responding width parameters have the following values 
P± ~ 5.3 fm and = 1-5 fm [23|. It should be pointed 
out that the THSR wave function Eq. ((2T|) is again of 
Gaussian type with a wide extension, centered at the 
origin. It is completely different from a Brink type wave 
function with the three a-particles placed at definite val- 
ues in space. A slight improvement of Eq. ^F2\ can still 
be achieved in taking the (ii parameters as Hill- Wheeler 
coordinates and superpose a couple of wave functions of 
the type ([22)) with different width parameters. Practi- 
cally 100% squared overlap with the wave function of the 
full RGM result of Kamimura et al. [11] is then achieved, 
as documented in [2^. It should be pointed out that 
the superposition of several Gaussians of the type ([2T|) 
does not at all change the physical content of the THSR 
wave function as a wide extended distribution centered 
around the origin. As we now have three a-particles, all 
in relative S'-states, one can begin to talk about coherent 
features, that is all a-particles occupying the same OS* 
orbit. Of course, the Pauli principle is acting, however 



weakly, and the OS* orbit is still occupied to over 70%, 
see Refs. [l^, [HI [111 ^^^^ discussion in Sec. VII. 

We also would like to attract the attention of the reader 
to the following important point: in spite of the fact that 
the THSR wave function describes ^Be very well, it is, 
of course, clear that no a-particle condensate aspect can 
be present with only two a-particles. This is also born 
out by the fact that in ®Bc the a-particle wave function 
still features quite strong influence of the Pauli principle 
with the two nodes seen on Fig.[T] On the other hand, as 
seen below on Fig. [H in Sect, lylll in and 1*^0, the 
a-particle wave functions in the condensate states have 
almost pure 05* wave character and, thus, the influence of 
the Pauli principle is much reduced (see Sec. HI) and the 
bosonic condensate feature born out. This stems from 
the fact that e.g. in the Hoyle state the a-particlcs arc, 
on average, by about 70 % (see Fig.O farther apart than 
in ^Be and also that the Hoyle state has to be orthogonal 
to the ground state of whereas the a-structure in *Be 
represents the ground state itself. 

Increasing the number of a-particles, the full RGM so- 
lution is not possible any longer. However, due to the rel- 
ative simplicity of the THSR ansatz, analogs to the Hoyle 
state have been found in ^^O, ^ONe, always situated close 
to the na disintegration threshold [l,[il] • Due to the high 
agreement with the full RGM results in the *Be and 
cases, one can expect that the THSR wave function also 
gives accurate results for the heavier systems, grasping 
well the physical situation of loosely bound a-particles 
moving in identical 05 orbits. Naturally the condensate 
aspect is realized the more, the larger the number of a's. 

Concluding this section one may just repeat the well 
known knowledge that the weakly bound na-particle 
states around the na disintegration threshold are not at 
all correctly described by standard a-cluster wave func- 
tions, with a crystal like structure of the a's. Rather 
the condensate aspect is dominant and imposes itself as 
the correct interpretation. For example the Hoyle state, 
therefore, can be seen as three almost inert a-particles 
moving in their own mean field potential, to good approx- 
imation given by a wide harmonic oscillator, whereas the 
a's are represented by four nucleons captured in narrow 
harmonic potentials. The situation is given as a cartoon 
in Fig. El 



III. INFLUENCE OF ANTISYMMETRISATION 

As already pointed out in the previous sections, a cru- 
cial question is whether for the Hoyle state the THSR 
wave function ((T|) together with (|4]) can be considered 
to good approximation as a product state of a particles 
condensed with their c.o.m. motion into the 05 orbital. 
For this, one has to quantify the influence of the anti- 
symmetriser ^ in ([1]). To a certain extent this question 
was already answered in earlier works [ll|. For instance 
in Ref. [ll| the single a particle density matrix, /o(7?., 7?.') 
was constructed and diagonalised for the Hoyle state. 
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FIG. 6: Pictorial representation of the THSR wave function 
for n = 3 (^■^C). The three a-particles are trapped in the 05*- 
state of a wide harmonic oscillator (B) and the four nucleons 
of each a are confined in the Os-state of a narrow one (b). All 
nucleons are antisymmetrised. 



The result was that the a's occupy the OS' orbit to over 70 
percent. The same result was later obtained with the so- 
called OCM (Orthogonality Condition Model) which is 
a very well tested method to describe cluster states [l^l ■ 
Because of the importance of this result in the present 
context, we present it again in Fig. [71 



1.0 



value of the antisymmetriser A, 

{B\A\B) 



{B\B) 



(23) 



where \B) is the THSR wave function ^ without the 
antisymmetriser, i.e. \'il)aiipa2'4'ai) see Eq. ([8]). The nor- 
malisation of the antisymmetriser A is chosen [3l| so that 
N{B) becomes unity in the limit where the inter-cluster 
overlap disappears, i.e. for the width parameter B ^ oo. 
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FIG. 7: Occupation of the single-Q orbitals of the Hoyle state 
of ^^C compared with the ground state [lol ]. 



We see that for the ground state the occupation num- 
ber distribution agrees with the SU{2i) shell model pic- 
ture (see Ref. [13 )> whereas for the Hoyle state the oc- 
cupation of the OS c.o.m. wave function of the a particle 
is, as mentioned, over 70 percent. It also is important 
to notice that no other state is occupied with more than 
about 7 percent meaning that the occupation of all other 
states is down by at least a factor of ten. This is a typical 
scenario for Bose condensed states. 

A more direct way to measure the influence of anti- 
symmetrisation is to consider the following expectation 



FIG. 8: The expectation value of the antisymmetriser for the 
product state \B). The values at the optimal B values, Bg 
for the ground and Bh for the Hoyle states, are denoted by 
circle and cross, respectively. See the text for the definition 
of \B), Bg and Bh- 

In Fig. [51 we show the expectation value N{B), 
Eq. of ^ as a function of the width parameter 

B. We chose, as optimal values of B for describing 
the ground and Hoyle states, B = Bg = 2.5 fm and 
B = Bh = 6.8 fm, for which the THSR states best 
approximate the ground state |0^) and the Hoyle state 
jOj'), respectively, which are obtained by solving the Hill- 
Wheeler equation. In fact, the normalised THSR state, 
|THSR(S))/^(THSR(B)|THSR(B)), gives the largest 
squared overlap 0.93 with the ground state |0]'') at i? = 
Bg. Similarly, it gives the largest squared overlap 0.78 
with the Hoyle state |0j) at i? = Bh- 

We should mention that Bg ^ b, since the ground 
state contains a-like correlations which lower the energy 
with respect to the limit of a pure Slater determinant 
{B = b = 1.35 fm) by roughly 5 percent [H, H^j- We see 
from Fig. [51 that N{Bh) ~ 0.62 and N{Bg) ~ 0.007 
what indicates that the influence of antisymmetrisation 
is strongly reduced in the Hoyle state compared with the 
one in the ground state. 

It is worth to have a closer look to the behavior of 
N{B)- It is seen that first this function raises very 
steeply, whereas for B > Bh the raise is much slowed 
down, reflecting the fact that the contribution from the 
one-nuclcon exchange term very slowly fades out [s^ . 
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FIG. 9: The squared overlap of the state \B) with the ground 
state and \Bg). The values at the optimal B values, Bg and 
Bh for the ground and Hoyle states, are marked by the circle 
and cross, respectively. See the text for the definition of \B), 
\Bg), Bg and Bh- 



An important point in the present considerations is 
that the THSR wave function for B ~ Bh is not auto- 
matically orthogonal to the ground state. This is con- 
trary to the situation with condensed cold bosonic atoms 
for which the density is so low that the overlap of the 
electron clouds can on average totally be neglected. It is 
nevertheless interesting to calculate the following overlap 
of I THSR) with the ground state |0^), obtained by solv- 
ing the Hill- Wheeler equation, or with \Bg), as a function 
of B, 

0^0+ 5. ^ l(0^|THSR)P 

^ 1 ' ' (THSR|THSR)' 
OiBg,B) 

|(THSR(gg)|THSR(B)}p 

~ (THSR(Bg)|THSR(Bg))(THSR(B)|THSR(B))' 

(24) 

From Fig. [HI we find that for both cases, the overlap is 
less than 0.12 indicating that orthogonality with ground 
state is nearly realised. 

In Fig. [ini the energy curves for the THSR wave func- 
tion with and without the antisymmetriser are shown, 
where 



E{B)^ 



(THSR|i?|THSR) _ {B\HA\B) 



(THSR|THSR) 
B\H\B) 



{B\A\B) 



{B\B) 



(25) 



They are measured from the ia threshold energy, E'tii 
?>Ea = —82.04 MeV, where is the binding energy 
of the intrinsic a particle [H, [2^, obtained with the 
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FIG. 10: The binding energies for |THSR(B)) and \B) as a 
function of B, denoted by E{B) and i5boson(B), respectively. 
See the text for the definition of E{B) and i5boson(B). They 
are measured from the calculated 3a threshold energy, Eth = 
3Ea = —82.04 MeV, where Ea is the binding energy of the 
intrinsic a particle. Volkov force No. 2 [s^l with Majorana 
parameter M = 0.59 and b = 1.35 fm are adopted as used in 
Refs. [2^, [2^. The values at the optimal B values, Bg and 
Bh for the ground and Hoyle states, are marked by the circle 
and cross, respectively. 



use of Volkov force No. 2 [33[. The second equality for 
E{B) in Eq. ([25]) holds due to the relation, [H,A] = 0. 
The minimum for E{B) is given ai B = Bg, which cor- 
responds to the ground state. The minimum energy 



E{Bg) - E, 



th 



-5.64 MeV, as also shown in Ref. 



On the other hand, i?boson(-B) — Eth ^ —100 MeV, gives 
unphysically large binding at small B values around Bg, 
indicating that antisymmetrisation plays an important 
role for the ground state. As B increases, however, the 
energy drastically gets smaller, and for B = Bh we have 
£^boson(^) — Eth = —17.5 MeV. This means that, com- 
pared with the ground state at B ~ Bg, the effect of 
antisymmetrisation is much reduced for the Hoyle state. 
It is, however, still essential to get the energy back on 
the spot. 

It is very important to point out that in this energy 
curve E{B) the second minimum corresponding to the 
Hoyle state is not present. This is because the THSR 
state with B = Bh, |THSR(Sjf)) stiU includes the 
ground-state components of about 10 percent what we 
have seen in Fig.[9l In fact, if we calculate the following 
binding energy. 



(P' 



(g.s. ) XHSRI H I P|*^-"- ^ THSR) 



THSR|P|^'-^THSR) 



(26) 



where the explicit orthogonalisation to the ground state 
is taken into account for the THSR state, with pj^ '' ' — 
1 — |0]'")(0j'"| like in Eq. (P^ . there appears the mini- 
mum corresponding to the Hoyle state at i? ~ Bh, as 
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FIG. 11: The binding energy in the orthogonal space to the 
ground state, denoted by Ep{B), together with E{B) in the 
previous figure. See the text for the definition of Ep{B). 
They are measured from the calculated 3a threshold energy, 
-Eth = -82.04 MeV. Volkov force No. 2 Q with Majorana 
parameter M = 0.59 and b — 1.35 fm are adopted as used in 
Refs. [2^, [2^. The values at the optimal B values, Bg and 
Bh for the ground and Hoyle states, are marked by the circle 
and cross, respectively. 



cnce of antisymmetrisation is expected in heavier Hoyle- 
like states. 



IV. DE BROGLIE WAVE LENGTH OF 
a-PARTICLES IN THE HOYLE STATE 




k [fm"'] 

FIG. 12: Momentum distribution of the a particle in the 
Hoyle state [lol ]. 



shown in Fig.[TTJ This is also discussed in Ref. [l^ where 
Brink-type 3a wave functions are used. We should also 
mention that the Hoyle state is much better approxi- 
mated by |P|S '-^THSR(S)) than by |THSR(B = Bh)), 
with proper normalisation factors for both states. The 
former state gives the largest squared overlap with the 
Hoyle state, 0.91 for B = 6.1 fm, which should be com- 
pared with 0.78 for the latter, a value already mentioned 
above. Thus the small admixture of the ground-state 
components for the THSR state is never negligible and 
the explicit elimination by pj^ '^ ' plays an essential role 
to describe the Hoyle state. It is thus true that the effect 
of the antisymmetrisation is not negligible even for the 
Hoyle state in a sense that the projection operator pj^ '' -' 
includes the compact ground-state components which are 
strongly subject to the antisymmetriser. Nevertheless, 
it is worth to emphasise that as a result of the explicit 
orthogonalisation to the ground state, the Hoyle state 
cannnot have a compact structure but has a dilute den- 
sity, for which, in the end, the effect of antisymmetrisa- 
tion is small. 

Therefore, let us point out again that the Hoyle state is 
to good approximation in a product state of three a par- 
ticles: about 70 percent of the a particles are in the 05 
orbit. Only less than 30 percent of the three alpha par- 
ticles are not in the bosonic product state, a fact due to 
antisymmetrisation. According to Ref. , it is expected 
that in heavier self conjugate nuclei the alpha particles 
arc still further apart because the stronger Coulomb re- 
pulsion lowers the Coulomb barrier. Thus even less influ- 



In this section we show, based on results of detailed 
microscopic calculations, that the de Broglie wave length 
A is larger by almost an order of magnitude than the inter 
a particle distance of about 3-4 fm. 

The de Broglie wave length of the a's moving in the 
Hoyle state can be estimated from the resonance energy 
of ^Be being roughly 100 keV. Else one can estimate the 
kinetic energy of the a-particles from a bosonic mean field 
picture using the Gross-Pitaevskii equation [sl] . Figure 3 
in Ref. [l^j shows the mean field potential of a-particles 
in the Hoyle state, with an indication of the position of 
the single a particle energy (180 keV). The kinetic en- 
ergy of the single a particle is calculated to be 380 keV. 
From this, the de Broglie wave length A = h/{2MaEy/^ 
is, therefore, estimated to be of a lower limit of approxi- 
mately 20 fm. A more reliable estimate of the de Broglie 
wave length is to use the expectation value of for the 
wave number k of the a particle in the Hoyle state, eval- 
uated from the momentum distribution of the alpha par- 
ticle, p{k), in Fig. [T^l obtained by a 3a OCM calcula- 
tion The resuh is A = 27r/y(fc2) ~ 20 fm, consis- 
tent with the previous value. These estimates all indicate 
that the de Broglie wave length is much longer than the 
inter a-particle distance, favoring a mean field approach. 

All these facts make the Bose aspects of the a's in 
the sense defined in the introduction, plausible and they 
may reveal specific features of coherence as implied by 
the notion of a bosonic product state. In the case of nu- 
clear pairing where the number of bosonic constituents, 
namely the Cooper pairs, is finite and not much larger 
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than the number of a's in light self conjugate nuclei, the 
collective properties of the 0+ ground states and the co- 
herence of these states have been revealed experimen- 
tally very early. The most conspicuous example being 
the strong reduction of the moment of inertia of super- 
fluid deformed nuclei from its classical value and the 
strong enhancement of the two-neutron transfer to the 
ground states 0, HI]. The even-odd staggering in the 
nuclear masses reveals pairing but not necessarily co- 
herence properties of the pairing state. All these effects 
of superfluidity are difficult to put into evidence for a- 
particle condensates for the simple reason that they are 
resonances around the a disintegration threshold with a 
finite life time. However, instead of pair-transfer, one 
can observe a-decay. Similar to two-neutron transfer, 
because of phase coherence, once a first alpha particle 
leaves the nuclear system, the probability that a second, 
third, ... will be emitted should be enhanced. This is 
precisely what we want to report in the following. 



V. MEASURABLE CONSEQUENCES OF 
LOOSELY BOUND ALPHA-PARTICLE STATES. 

It is of great importance and interest to discuss even- 
tual measurable consequences for the signature of boson- 
condensates in nuclei. So far practically all measured 
quantities of the Hoyle state in ^^C have been repro- 
duced with the THSR wave function with rather good 
accuracy without any adjustable parameter poj . For in- 
stance it can be affirmed that the Hoyle state has quite 
dilute density, about only 1/3 ^ 1/4 of the one of the 
ground state of ^^C. The density of the Hoyle state is 
about the same, or slightly less than the one of *Be. To 
illustrate this dramatic effect it is sufficient to state that 
the density at the origin is only half of the one in the 
ground state of ^^C This fact alone suggests that 

in the Hoyle state and other similar states in heavier na 
nuclei mentioned above, there is enough space that nucle- 
ons cluster into three alphas. The latter mostly interact 
via the *Be resonance. 

We now will give new interpretations of older exper- 
iments involving multiple a-particle decay out of com- 
pound states which reveal coherence effects of an na- 
particle gas. The idea is that in excited states of heavier 
N = Z compound nuclei above a certain excitation en- 
ergy a low density state of a-particles can be formed, 
in a mixed phase of fermions and bosons [36| . Such 
states can also be formed on top of an inert core like, 
e.g. ^^O or ^''Ca. The decay process is expected, in the 
light of the present considerations, to show special fea- 
tures as exemplified below. In particular the multiple 
Qf-decay observed does not correspond to results of the 
Hauser-Feshbach theory of compound nuclear decay. 

Previous studies of ^Be-emission from excited com- 
pound nuclei [13, IH, with particle-7-coincidcnces 
have shown strong effects in the 7-spectra, if statistical a- 
a emission was compared with the *Be emission. In these 



experiments the multiple emission of a-particles is reg- 
istered with the ISIS-Ai?-i? particle detection array, in 
coincidence with the 7-decay of the residual nuclei regis- 
tered with the 7-detector array GASP at the Laboratorii 
Nationale di Legnaro (Italy). As an example, we cite the 
reaction ^SSi + ^^Mg ^^Fc at 130 MeV ^ forming a 
compound nucleus at the excitation energy of = 76 
MeV. For the a-a-corrclations in these experiments an 
enhancement is observed when they are emitted into the 
same direction. In view of the opening angle (27 degr.) 
of one individual telescope of the ISIS-Ai?-i?-system (42 
telescopes), we were able to detect the spontaneous de- 
cay of the unbound states just at the decay-thresholds, 
namely of the ®Be and the "'^^C*(02') states, into two or 
three a-particles, respectively. The decay energies give 
relative energies of less than 100 keV and opening angles 
result, which fit into the solid angle of a single telescope. 
The corresponding coincident (particle gated) 7-decays 
are compared with the spectra obtained from statistical 
a-particle emission into different telescopes, but with the 
same a-multiplicity. 

This particularly striking effect is illustrated in Fig. [T3] 
in the case of the above mentioned reaction ^®Si-|-^^Mg — > 
^^Fe — > ^°Ca-|-3a. We compare the -'^^C*-emission (lower 
spectrum shown in Fig. I13p . with that triggered by the 
statistical decay with three random a-particles (here the 
7 spectra are dominated by transitions in ^°Ca and ■^^K). 
Quite conspicuous additional gamma rays of '^^Ar appear 
in the spectrum gated with ^^C* -emission, implying that 
a 4th a-particle is emitted, which is not predicted in the 
Hauser-Feshbach approach for statistical compound de- 
cay. 

Further explanation of the observed effect within the 
concept of a gas of almost ideal bosons in ^^Fe, (in this 
case with a ■**'Ca-core), has been proposed in Ref. [40| . 
If such a compound state is formed in the cited reaction, 
its characteristic feature is a large radial extension with 
a very large diffuseness of the density distribution. With 
the large diffuseness of the potential the calculation for 
the emission of the ^^C* resonance gives a dramatic low- 
ering of the emission barrier, as compared to the statis- 
tical multi-a-particle emission, it amounts to more than 
10 MeV. Thus the residual nucleus (^°Ca) is populated 
at much higher excitation energy, and a further a-dccay 
can occur. 

The coherent properties of the threshold states consist- 
ing of a-particles interacting via their resonances (36j , are 
due to properties of the ®Bc ground state or the ^^C*(Oj)- 
state. As reported in Sect. llVi the a-particles have a 
large dc Broglic wave length of relative motion, of A = 
/i/(2M„£;)i/2 _ 20 fm, or more. In the decay of the 
compound nucleus the decay steps are usually statisti- 
cally independent and, thus, the a particles leave the 
nucleus one by one. However, for the coherent state the 
a's are already existing (large spectroscopic factors) and 
their wave functions overlap coherently, the decay may 
be a simultaneous na decay, keeping the phase relations 
of their relative motion. In this way also the Coulomb 



11 




2400 

E [keV] 



FIG. 13: Coincident 7-spectra gated with Ai?-_E-signals with the emission of three random q's in different detectors, in 
comparison with that obtained by the ^'^C* (Oj)-gate, as indicated [s^. The reaction is ^**Si+^''Mg at 130 MeV. Note the 
additional '^^Ar-line in the lower spectrum, which indicates an additional emission of one a. Figure is adopted from Ref. [38l |. 



barrier gets lowered, a fact which further enhances the 
decay probability. This leads to the observed unbound 
resonances of ^Be and ^^C*(Oj). This result can be in- 
terpreted as the observation of bosonic coherence in the 
compound nucleus. 

Such enhanced decays also have been observed in stud- 
ies of ^Be-emissionjSa, [1^ SJ] with the corresponding 
gated 7-spectra. When compared with the statistical de- 
cays into two alpha-particles additional gamma transi- 
tions are observed. At the time of these experiments, 
the discussion of a-particle coherence was not consid- 
ered, attempts to explain the observed effects within the 
extended Hauser-Feshbach (EHF) formalism failed [3^. 

To summarize we can state that future dedicated ex- 
periments, with Ai^-iJ-telescopes in coincidence with 
an efficient 7-detection array as described, may be well 
suited to establish the existence of THSR type states in 
excited N = Z nuclei. 



VI. DECAY PROPERTIES 

One may ask the question whether Hoylc-likc states 
in nuclei heavier than ^^C can exist. An argument can 
be based on the fact that the alpha-particle conden- 
sate states occur near the alpha-particle disintegration 
threshold which rapidly grows in energy with mass and 
thus the level density in which such a condensate state 
is embedded raises enormously. For example the alpha- 
disintegration threshold in ^^C is at Ex = 7.27 MeV and 
in ^^O it is already at ~ 14.4 MeV. Under ordinary 
circumstances this could mean that the alpha-particle 
THSR state in ^^O, which we suppose to be the well 
known 0+ state at E^ = 15.1 MeV [l3l, has a very short 
life time and in Ref. a Fermi gas estimate is made 
in this respect. However, on the one hand it is a fact 
that the supposed ^^O "Hoyle" -state at E^ ^ 15.1 McV 
has experimentally, for such a high excitation energy, a 
startling small width of 160 keV and on the other hand 
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it is easily understandable that such an exotic configura- 
tion as four alpha-particles moving almost independently 
within the common Coulomb barrier, has great difficul- 
ties to decay into states lower in energy which all have 
very different configurations. How else one could explain 
such a small width of a state this high up in energy? 
It is precisely one of the strong indications of Hoyle-likc 
states that they should be unusually long lived! It is 
furthermore well known that the Hoyle state cannot be 
explained even with the most advanced shell model cal- 
culations. Its energy comes at 2-3 times its experimen- 
tal value [l^l . This is a clear indication that shell model 
configurations only couple extremely weakly to alpha gas 
states. One can argue that many of the states in ^^O be- 
low Ex = 15.1 MeV are of shell model type. There are 
also alpha- ^^C configurations but since ^^C also has shell 
model configuration, it again is difficult for the four alpha 
condensate state to decay into. 



decay- 
ing into possible channels and the total width. The reduced 
widths defined in Eq. (|28|) are also shown, a is the channel 
radius. 



[a = 8.0 fm) (g = 7.4 fm) (a = 8.0 fm) 



Tl (keV) 



26 
0.006 



0.004 



2 X 10" 
0.15 



Total 



34 



Let us make a more quantitative estimate of the decay 
width of the Ej.. = 15.1 MeV state. Based on the R- 
matrix theory [43|, the decay width Tl can be given by 
the following fomulae: 



ka 

^^^"^ = Fl{ka) + Gl{kay 
7i(a)=0i(a)7^(a), 

7w(a) 



2/ia2 



(27) 



where k, a and fi are the wave number of the relative 
motion, the channel radius, and the reduced mass, re- 
spectively, and Fl, Gl, and Pl{o-) ai'c the regular and 
irregular Coulomb wave functions and the corresponding 
penetration factor, respectively. The reduced width of 
Oj^{a) is related with the wave function '5(0^) of the al- 
pha condensate in ^^O obtained as the sixth 0+ state in 
Ref. [lal, as follows: 



5{r' - a) 



YL{r')^LrC) ^ vl/(0+) 



(28) 



where ^^(^^C) is the wave function of ^^C, given by the 
3a OCM calculation [l^. In Table [H we show the par- 
tial a decay widths of the Ot state Tl decaying into the 



a-hi2C(0+), a-hi2C(2+) and a-h^^ c(0+) channels, total 
a decay width which is obtained as a sum of the partial 
widths, and reduced widths 0'j^{a) defined in Eq. (|28|1 . 
Experimental values are all taken as given by the decay 
energies. Thus the excitation energy of the calculated 0^ 
state is assumed to be 15.1 MeV, the one of the observed 
0,^ state. 

The obtained very small total a decay width of 34 keV, 
in reasonable agreement with the corresponding experi- 
mental value of 160 keV indicates that this state is unusu- 
ally long lived. The reason of this fact can be explained in 
terms of the present analysis as follows: Since this state 
has very exotic structure composed of gas-like four alpha 
particles, the overlap between this state and a+^^ C{0^) 
or a C{2^) wave functions with a certain channel 
radius becomes very small, as this is, indeed, indicated 
by small Oj^{a) values, 0.006 and 0.004, respectively, and 
therefore by small 7£(a) values. These largely suppress 
the decay widths expressed by Eq. ([27|) in spite of large 
values of penetration factors caused by large decay en- 
ergies 7.9 MeV and 3.5 MeV into these two channels, 
a C(Of ) and a-h^^ C(2:^), respectively On the other 
hand, the decay into a +^'^ C(Oj) is also suppressed due 
to very small penetration caused by very small decay 
energy 0.28 MeV into this channel, even though the cor- 
responding reduced width takes a relatively large value 
02 (a) = 0.15, which is natural since the Oj state of ^^c 
has a gas-like three-alpha-particle structure. It is very 
likely that the above mechanism holds generally for the 
alpha gas states in heavier na systems, and therefore such 
states can also be expected to exist in heavier systems as 
a relatively long lived resonances. 



VII. SIMILARITY OF a-PARTICLE WAVE 
FUNCTIONS IN HOYLE-LIKE STATES 



In Fig. [T3] we show, side by side, radial parts of the 
single-a S orbits (for definition, see Refs. [13, [Hjlill) of 
the Hoyle state (^^c) and the 0^ state in ^^O. We see 
an almost identical shape. Of course, the extension is 
slightly different because of the smallness of the system, 
i.e., we are not dealing with a macroscopic condensate as 
is discussed above. The nodeless character of the wave 
function is very pronounced and only some oscillations 
with small amplitude arc present in ^^C, reficcting the 
weak infiuence of the Pauli principle between the a's, see 
the discussion in Sec. II. On the contrary, due to their 
much reduced radii, the "a-like" clusters strongly over- 
lap in the ground states of ^^C and ^^O, producing strong 
amplitude oscillations which take care of antisymmetri- 
sation between clusters [l3, 113 • This example demon- 
strates the bosonic product nature of the Hoyle state and 
the state in ^^O. 
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FIG. 14: Radial distributions of the single-a S orbits, (a) of the state in ^^C (Hoyle state) and (b) of the state in ^^O. 



VIII. DISCUSSIONS, SUMMARY, AND 
CONCLUSIONS 

In this work we considerably deepened several aspects 
of the THSR description of low density a-particle states 
in self-conjugate nuclei. We show that the THSR wave 
function which has alpha-particle condensate structure 
in analogy with the number projected BCS wave func- 
tion for Cooper pairs, grasps the physics of loosely bound 
alpha-particle states as the ^Bc ground state and the 
Hoyle state in ^^C, etc. much better than the usual 
cluster wave functions of the Brink type where rather a 
crystal structure is involved, the alpha particles with free 
space extension being placed at certain geometrical posi- 
tions with respect to one another in the nucleus. Indeed, 
wc have shown for the example of ^Bc that the superpo- 
sition of about 30 Brink type wave functions is needed to 
describe ®Be ground state with the same accuracy as the 
single component THSR wave function which practically 
coincides with the exact solution of the RGM wave func- 
tion. Similar results are obtained concerning the Hoyle 
state, where about 55 components of the Brink type are 
needed 

One entire section is dedicated to the study of effects 
from antisymmetrisation between the alpha particles in 
the THSR (Hoyle) state. We studied the expectation 
value N{B) of the antisymmetriscr as a function of the 
width parameter B in the THSR wave function which de- 
termines the mean distance between the a-particles and 
found that N{B) raises very fast as a function of B. For 
the Hoyle state with B = Bh, N{Bh) is about 0.62 
whereas for the ground state, i.e. B = Bg we find N{Bg) 
is about 0.007. So this value increases from the ground 
state to the Hoyle state by about a factor of hundred, 
indicating the strongly reduced action of antisymmetri- 
sation in the Hoyle state. Similar conclusions are found 
for energies and orthonormality relations, sec Figs. ISlfTT] 

We also discussed the possible experimental conse- 
quences of alpha-particle coherence. This question is 
more difficult to answer than for nuclear pairing by the 



fact that those states appear only at low density {p = 
Po/3 ~ Po/4) and thus correspond to excited states be- 
ing of unusual long life time on nuclear scales 10~^^ s 
for e.g. the Hoyle state). Transfer experiments of alpha- 
particles or measurements of moments of inertia which so 
clearly demonstrate superfluid features of nuclei in the 
case of pairing are, therefore, very difficult to conceive 
in the case of alpha-particles. However, instead of trans- 
fer, one may investigate decay. For example in an excited 
state of ^^Fe, if there exists coherence of a certain number 
of alpha-particles on top of an inert core, the simultane- 
ous decay of two or more alpha-particles will be enhanced 
with respect to a purely statistical decay. Exactly this 
feature has been observed and related to alpha-particle 
condensation in the reaction ^^Si -|- ^"'Mg at E'lab = 130 
McV. However, other experiments may be conceived in 
the future. For example, with a heavy ion reaction e.g. 
^*Si may be (Coulomb) excited to the Ikeda threshold 
of seven alpha break up, that is = 38.46 MeV [i^ 
and then seven alpha-particles may expand as a coher- 
ent state verifiable with performant multi-particle detec- 
tors (46j . Measuring energies and angles of the a's may 
allow to establish an invariant mass spectrum which iden- 
tifies THSR states, even if they are relatively broad and 
hidden among other states, seen via gammas or particle 
evaporation. 

Since the alpha-particle condensate states appear 
around the na decay threshold, these states are the 
higher up in the continuum, the heavier the nuclei be- 
come. For example in ^^C the Ikeda threshold for al- 
pha decay [H] is at = 7.27 MeV, in ^''O it is at 
E^ = 14.4 MeV, in ^^Mg at E^ = 28.48 McV and so 
on. One may object that, because of their high excita- 
tion energies, those analogs to the Hoyle state will decay 
very quickly. However, we argued that these alpha gas 
states have very unusual structure and thus couple only 
very weakly to states at lower energy. The example of 
the Hoyle state shows that all states of shell model struc- 
ture are practically decoupled. This makes up for the 
large majority of states. For the supposedly analog to 
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the Hoyle state, namely the 0^ state at = 15.1 MeV 
in ^^O, a quantitative estimate of the decay was made, 
explaining the very small width of 160 keV. 

Going on with this reasoning, it is not at all excluded 
that alpha gas states in even heavier self conjugate nu- 
clei will have a quite unusual long life time, given their 
excitation energy high up in the continuum. 

Furthermore, the condensate character of the alpha- 
gas states also has been pointed out in showing, see 
Fig. [131 that the condensate wave function of one alpha- 
particle changes, apart from a trivial size effect, very lit- 
tle in going e.g. from the Hoyle state in ^^C to the cor- 
responding state at — 15.1 MeV in ^^O, that is, a 
condensate wave function is a OS" state. 

One of our strong arguments that the alpha-particles 
in the Hoyle state and possibly in the 0^ of ^^O, form 
an alpha-particle gas, captured inside the Coulomb bar- 
rier, is deduced from the fact that we constructed a sin- 
gle alpha-density matrix whose eigenvalues show that the 
alpha product states are realised to around 60-70% , al l 
bosons being in the lowest quantum state [l3, [HI, Il3 |. 
The occupancies of all higher quantum levels are down by 
at least a factor of ten. However, the authors of Ref. [3 
mentioned that the way how to define the density matrix 
in a self bound Bose system with a finite number of par- 
ticles, is not unique and that different definitions might 
give different occupancies. We followed in this respect 
the line of thought outlined by Pethick and Pitaevskii 
in Refs. [13, SI] where they say that if in a homogeneous 
system there is Bose condensation, then there is no reason 
that, if the same system is put into an external potential 
or if the system is self bound in a mean field potential, the 
system be not also in a non-fragmented condensate state, 
as long as the intrinsic system is not excited. We showed 
that our definition of the boson density matrix satis- 
fies this physically very reasonable boundary condition 
in using Jacobi coordinates for the internal system [i^ . 
We furthermore showed that for the OS* harmonic oscil- 
lator wave function the internal one-body density matrix 
is uniquely determined under another reasonable condi- 
tion [49| . The uniqueness for more general wave functions 
is also demonstrated [13, [Ml . 

In the light of this finding, we would like to discuss 
again the content of the THSR alpha-particle condensate 
wave function Eq. It is very important to remark, as 
is explained in Ref. [1[ , that this antisymmetrised a parti- 
cle product wave function contains two limits exactly. On 
the one hand, for B = b we have a pure harmonic oscilla- 
tor wave function because the antisymmetriser generates 
out of the product of simple Gaussians all hig her nodal 
wave functions of the harmonic oscillator |52l |. On the 
other hand, for B ^ b the THSR wave function tends to 
a pure product state of alpha-particles, i.e. a mean field 
wave function, since in this case the antisymmetriser can 
be neglected. Indeed B triggers the extension of the nu- 



cleus, i.e. its average density. For alpha particles kept 
at their free space size (small 6), the alpha-particles are 
then for large S-values far apart from one another and 
do not feel any action from the Pauli principle, see a de- 
tailed discussion of the action of the antisymmetriser as 
a function of density in Ref. [10]. The question is then 
whether, e.g. for the Hoyle state, the above wave func- 
tion is closer to a shell model like Slater determinant or to 
an a-particle product state. Precisely this question is an- 
swered by the above discussed eigenvalues of the density 
matrix. In this respect it is important to point out that 
in the calculation of the afore-mentioned density matrix 
always the total c.o.m. motion has been split off in the 
wave function of Eq. ([T]) and that for the remaining rel- 
ative c.o.m. coordinates the Jacobi ones have been used, 
as is clearly explained in Refs. [l3j[li|- In Refs. (Tol. [Tl| 
it has been shown, as explained, that the a's in the Hoyle 
state occupy to over 70 % the OS'-orbit. Therefore, the 
Hoyle state is in good approximation a product of three 
alpha particles, that is a condensate. This finding also is 
corroborated by our study on antisymmetrisation effects 
in the Hoyle state in Sec. HI. As already mentioned, it 
has been found that the effects of antisymmetrisation are 
weak. 

In summary we can say that our study clearly shows 
that the loosely bound alpha-particle states of very low 
density, close to the decay threshold in self-conjugate nu- 
clei, are characterised by a shallow self consistent mean 
field of wide extension in which the c.o.m. motion of the 
alpha-particles occupies with about 70% the lowest OS 
level. In spite of the very different number of particles 
and other important differences, the situation has, there- 
fore, some analogy with the case of cold atoms. Avoid- 
ing vague and qualitative arguments, we hope to have 
been sufficiently detailed and convincing, based on new 
insights from precise numerical results, to confirm the 
existence of low density bosonic a particle gas states in 
nuclei and to affirm the usefuUness of this novel concept 
in nuclear physics described with the THSR wave func- 
tion. 



Acknowledgements 

The authors express special thanks to T. Neff whose 
detailed studies of the ®Be wave function, made avail- 
able to us, helped, together with extended discussions, 
our work appreciably. They also thank Kanto Gakuin 
University (KGU) and Yukawa Institute for Theoretical 
Physics at Kyoto University, Japan. Discussions during 
the KGU Yokohama Autumn School of Nuclear Physics 
and the YIPQS international molecule workshop held in 
October 2008 were useful to complete this work. 



15 



[1] A. Tohsaki, H. Horiuchi, P. Schuck and G. Ropke, Phys. 
Rev. Lett. 87, 192501 (2001). 

[2] N.T. Zinner and A.S. Jensen, Phys. Rev. C 78, 
041306(R) (2008). 

[3] G. Ropke, A. Schnell, P. Schuck, and P. Nozieres, Phys. 
Rev. Lett. 80, 3177 (1998). 

[4] M. Beyer, S.A. Sofianes, C. Kuhrts, G. Ropke, and P. 
Schuck, Phys. Lett. 448B, 247, (2000). 

[5] T. Sogo, R. Lazauskas, G. Ropke, and P. Schuck, Phys. 
Rev. C 79, 051301(R) (2009). 

[6] B. Doucot, J. Vidal, Phys. Rev. Lett. 88, 227005 (2002); 
S. Capponi, G. Roux, P. Lecheminant, P. Azaria, E. 
Boulat, S.R. White, Phys. Rev. A 77, 013624 (2008); 
P. Lecheminant, E. Boulat, and P. Azaria, Phys. Rev. 
Lett. 95, 240402 (2005). 

[7] A. Bohr and B. R.Mottelson, Nuclear Structure (Ben- 
jamin, New York, 1975), Vol. 2; P. Ring and P. Schuck, 
The Nuclear Many Body Problem (Springer- Verlag, New 
York, 1980); J. -P. Blaizot and G. Ripka, Quantum The- 
ory of Finite Systems (MIT, Cambridge, MA, 1986). 

[8] For example, P. Schuck, Y. Funaki, H. Horiuchi, G. 
Ropke, A. Thosaki, and T. Yamada, Prog. Part. Nucl. 
Phys. 59, 285 (2007). 

[9] A. Tohsaki, H. Horiuchi, P. Schuck and G. Ropke, Nucl. 
Phys. A 738, 259 (2004). 
[10] T. Yamada and P. Schuck, Eur. Phys. J. A. 26, 185 
(2005). 

[11] H. Matsumura and Y. Suzuki, Nucl. Phys. A 739, 238 
(2004). 

[12] Y. Funaki, T. Yamada, H. Horiuchi, G. Ropke, P. Schuck, 

and A. Tohsaki, Phys. Rev. Lett. 101, 082502 (2008). 
[13] D.M. Brink, in Proceedings of the International School 

of Physics "Enrico Fermi", Course 36 (Academic Press, 

New York, London, 1966) p. 247. 
[14] R. Tamagaki et ai, Prog. Theor. Phys. Suppl. 52, 1 

(1972), and references therein. 
[15] For example, Y. Fujiwara, H. Horiuchi, K. Ikeda, M. 

Kamimura, K. Kato, Y. Suzuki, and E. Uegaki, Prog. 

Theor. Phys. Suppl. 68, 29 (1980). 
[16] T. Yamaya, K. Katori, M. Fujiwara, S. Kato and S. 

Ohkubo, Prog. Theor. Phys. Suppl. 132, 73 (1998); F. 

Michel, S. Ohkubo and G. Reidemeister, Prog. Theor. 

Phys. Suppl. 132, 7 (1998). 
[17] T. Wada and H. Horiuchi, Phys. Rev. C 38, 2063 (1988). 
[18] N. Itagaki, T. Otsuka, K. Ikeda and S. Okabe, Phys. Rev. 

Lett. 92, 142501 (2004). 
[19] J. A. Wheeler, Phys. Rev. 52, 1083 (1937); ibid., 1107 

(1937) 

[20] K. Ikeda, R. Tamagaki, S. Saito, H. Horiuchi, A. Tohsaki- 
Suzuki and M. Kamimura, Prog. Theor. Phys. Suppl. 62 
(1977). 

[21] R. B. Wiringa, S. C. Pieper, J. Carlson, and V. R. Pand- 
haripande, Phys. Rev. C 62, 014001 (2000). 

[22] Y. Funaki, H. Horiuchi, A. Tohsaki, P. Schuck and G. 
Ropke, Prog. Theor. Phys. 108, 297 (2002). 

[23] The 2a wave function given in Ref. [22| as the exact 2a 
RGM solution is different from the one given here. We 
now realized that to a better solution, i.e. to a better 
energy convergence within a hundred eV, we should also 
include the components of Brink wave function with val- 
ues of R as far as 30 fm, while previously the Brink com- 



ponents only up to i? = 12.5 fm have been adopted. 

[24] Y. Funaki, H. Horiuchi, G. Ropke, P. Schuck, A. Tohsaki 
and T. Yamada, Phys. Rev. C 77, 064312 (2008). 

[25] E. Uegaki, S. Okabe, Y. Abe, and H. Tanaka, Prog. 
Theor. Phys. 57, 1262 (1977); E. Uegaki, Y. Abe, S. Ok- 
abe, and H. Tanaka, Prog. Theor. Phys. 59, 1031 (1978); 
62, 1621 (1979). 

[26] M. Chernykh, H. Feldmeier, T. Neff, P. von Neumann- 
Cosel, and A. Richter, Phys. Rev. Lett. 98, 032501 
(2007). 

[27] The reason why the present value for the squared overlap 
is different from 0.97 which was given in Ref. [29| is as 
follows: In Ref. [2^, as the ground state configuration 
used in the projection operator P^''''\ we took the wave 
function which corresponds to the minimum of the en- 
ergy surface in Fig. 1 of Ref. [2^ . However in the present 
definition of P'f'^'\ we use the ground-state solution of 
the Hill- Wheeler equation in the two-parameter space of 
and /3z. 

[28] Y. Fukushima and M. Kamimura, Proc. Int. Conf. on 
Nuclear Structure, Tokyo, 1977, ed. T. Marumori (Suppl. 
of J. Phys. Soc. Japan, 44, 225 (1978)); M. Kamimura, 
Nucl. Phys. A 351, 456 (1981). 

[29] Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck and G. 
Ropke, Phys. Rev. C 67, 051306(R) (2003). 

[30] S. Saito, Prog. Theor. Phys. 40, 893 (1968); 41, 705 
(1969); Prog. Theor. Phys. Suppl. No. 62, 11 (1977). 

[31] A. Tohsaki, Y. Funaki, H. Horiuchi, G. Ropke, P. Schuck 
and T. Yamada, in Proceedings of KGU Yokohama Au- 
tumn School of Nuclear Physics, Yokohama, 2008, ed. T. 
Yamada and Y. Funaki (Int. J. Mod. Phys. A 24, 2003 
(2009)). 

[32] A. Tohsaki-Suzuki, Prog. Theor. Phys. Suppl. 62, 191 
(1977). 

[33] A.B. Volkov, Nucl. Phys. 74, 33 (1965). 
[34] T. Yamada and P. Schuck, Phys. Rev. C 69, 024309 
(2004). 

[35] W. von Oertzen and A Vitturi, Rep. Prog, in Phys. 64, 
1247 (2001). 

[36] W. von Oertzen, Eur. Phys. J. A 29, 133 (2006). 

[37] W. von Oertzen, Phys. Scr. T 88, 83 (2000). 

[38] Tz. Kokalova et ai, Eur. Phys. J. A 23, 19 (2005). 

[39] S. Thummerer, W. von Oertzen, B. Gebauer, S. Lenzi, 

A. Gadea, D.R. Napoli, C. Beck and M. Rousseau, J. 

Phys. G 27 1405 (2001). 
[40] Tz. Kokalova, N. Itagaki, W. von Oertzen, and C. Whel- 

don, Phys. Rev. Lett. 96, 192502 (2006). 
[41] S. Thummerer et ai, Phys. Scr. T 88, 114 (2000). 
[42] B. R. Barrett, B. Mihaila, S. C. Pieper, and R. B. 

Wiringa, Nucl. Phys. News, 13, 17 (2003), and B. R. 

Barrett, private communication. 
[43] A. M. Lane and R. G. Thomas, Rev. Mod. Phys. 30, 257 

(1958). 

[44] T. Yamada, Y. Funaki, H. Horiuchi, G. Ropke, P. Schuck, 
and A. Tohsaki, Phys. Rev. A 78, 035603 (2008). 

[45] K. Ikeda, N. Takigawa, and H. Horiuchi, Prog. Theor. 
Phys. Suppl. extra number, 464 (1968). 

[46] I. Tanihata and W. von Oertzen, private communication. 

[47] C.J. Pethick and L.P. Pitaevsku, Phys. Rev. A 62, 
033609 (2000). 



16 



[48] L.P. Pitaevskii, private communication. 

[49] T. Yamada, Y. Funaki, H. Horiuchi, G. Ropke, P. Schuck, 

and A. Tohsaki, Phys. Rev. C 79, 054314 (2009). 
[50] Y. Suzuki and M. Takahashi, Phys. Rev. C 65, 064318 

(2002). 



[51] Y. Suzuki, W. Horiuchi, M. Orabi, K. Arai, Few-Body 

Systems 42, 33 (2008). 
[52] T. Yamada, Y. Funaki, H. Horiuchi, K. Ikeda, and A. 

Tohsaki, Prog. Theor. Phys. 120, 1139 (2008). 



